DESY05-110 
July 2005 



ISSN 0418-9833 



Measurement of F^^ and F^^ at Low and x 
using the HI Vertex Detector at HERA 



HI Collaboration 



Abstract 

Measurements are presented of inclusive charm and beauty cross sections in e~^p collisions 
at HERA for values of photon virtuality 12 < < 60 GeV^ and of the Bjorken scaling 
variable 0.0002 < x < 0.005. The fractions of events containing charm and beauty quarks 
are determined using a method based on the impact parameter, in the transverse plane, of 
tracks to the primary vertex, as measured by the HI vertex detector. Values for the structure 
functions F^^^ and are obtained. This is the first measurement of in this kinematic 
range. The results are found to be compatible with the predictions of perturbative quantum 
chromodynamics and with previous measurements of ^2*^. 
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1 Introduction 



Measurements of the charm (c) and beauty (b) contributions to the inclusive proton structure 
function F2 have been made recently in Deep Inelastic Scattering (DIS) at HERA, using infor- 
mation from the HI vertex detector, for values of the negative square of the four momentum of 
the exchanged boson Q'^ > 150 GeV^ [1]. In this high Q'^ region a fraction of ~ 18% (~ 3%) 
of DIS events contain c (6) quarks. It was found that perturbative QCD (pQCD) calculations at 
next-to-leading order (NLO) gave a good description of the data. In this paper a similar method 
is employed, using data from the same running period, to extend the measurements to the range 
of lower Q^, 12 < < 60 GeV^ and of Bjorken x, 0.000197 < a: < 0.005. 

Events containing heavy quarks are distinguished from those containing only light quarks 
by reconstructing the displacement of tracks from the primary vertex, using precise spatial 
information from the HI vertex detector. The long lifetimes of c and b flavoured hadrons lead 
to larger displacements than for light quark events. The charm structure function and the 
beauty structure function are obtained from the measured c and b cross sections after small 
corrections for the longitudinal structure functions and F^^. The measurements at low 
Q"^ benefit from increased statistics when compared to those at high Q'^. However, the low 

region is experimentally more challenging because the final state does not have as large a 
transverse boost in the laboratory frame. The separation between b and c events is also difficult 
since, although the c fraction is expected to be similar as at high Q^, the b fraction is expected 
to be much smaller (~ 0.6% at = 12 GeV^ [2,3]). 

Previous measurements of the open charm cross section in DIS at HERA have mainly been 
of exclusive D or D* meson production [4-6]. From the D* measurements the contribution of 
charm to the proton structure function has been derived by correcting for the fragmentation frac- 
tion /(c — i> D*) and the unmeasured phase space (mainly at low values of transverse momentum 
of the meson). The results are found to be in good agreement with pQCD predictions. The b 
cross section in DIS, in a similar kinematic region to the present analysis, has been measured 
for events containing a muon and an associated jet in the Breit frame in the final state [7, 8]. 
The measured cross sections are found to be somewhat higher than perturbative calculations at 
NLO. 

2 Theoretical Description of Heavy Flavour Production in 
Deep Inelastic Scattering 

2.1 NLO QCD Calculations 

In the framework of NLO QCD analyses of global inclusive and jet cross section measure- 
ments, the production of heavy flavours is described using the variable flavour number scheme 
(VFNS) which aims to provide reliable pQCD predictions over the whole kinematic range in 
Q"^. At values of ^ the effects of the quark mass m must be taken into account and the 
heavy flavour partons are treated as massive quarks. The dominant LO process in this region is 
photon gluon fusion (PGF) and the NLO diagrams are of order [9]. As Q'^ increases, in the 
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region ^ m^, the heavy quark may be treated as a massless parton in the proton. Several 
approaches [10-12] have been developed which deal with the transition from the heavy quark 
mass effects at low Q'^ to the asymptotic massless parton behaviour at high Q'^. Recently, pre- 
dictions for inclusive heavy flavour production within a VFNS approach have been calculated 
at next-to-next-to-leading order (NNLO) [13]. 

Predictions for the charm and beauty cross sections may also be obtained from fits [14] to 
the HERA inclusive F2 data based on CCFM evolution [15]. The heavy quarks are produced 
in the fixed flavour number scheme (FFNS) according to the LO PGF off-shell matrix elements 
(with rric = 1.5 GeV and = 4.75 GeV) convoluted with the CCFM A;r-unintegrated gluon 
density of the proton (J2003 set 1 [14]). The predictions are calculated using the Monte Carlo 
program CASCADE [16]. 

2.2 IMonte Carlo Simulation 

Monte Carlo simulations are used to correct for the effects of the finite detector resolution, 
acceptance and efficiency. The Monte Carlo program RAPGAP [17] is used to generate low 
DIS events for the processes ep — > ebbX and ep — > eccX. The Monte Carlo program 
DJANGO [18] is used to generate light quark (uds) events. Both programs combine 0(as) ma- 
trix elements with higher order QCD effects modelled by the emission of parton showers. The 
heavy flavour event samples are generated according to the massive PGF matrix element with 
the mass of the c and b quarks set to rric, = 1.5 GeV and rrih = 4.75 GeV, respectively. In the 
heavy flavour event generation, the DIS cross section is calculated using the parton distribution 
functions (PDFs) from [19]. The light flavour event samples are generated with the LO PDFs 
from [20]. The partonic system for all generated events is fragmented according to the LUND 
string model implemented within the JETSET program [21]. The HERACLES program [22] 
calculates single photon radiative emissions off the lepton line, virtual and electroweak correc- 
tions. The Monte Carlo program PHOJET [23] is used to simulate the background contribution 
from photoproduction (7^ — > X). 

The samples of events generated for the uds, c, and b processes are passed through a detailed 
simulation of the detector response based on the GEANT3 program [24], and through the same 
reconstruction software as is used for the data. A total of 50 million uds events, 9 million c 
events and 1 million b events were simulated to evaluate the cross sections, corresponding to 
luminosities of 90 pb"^, 160 pb"^ and 980 pb"^, respectively. 

3 HI Detector 

The analysis is based on a low sample of e~^p neutral current scattering events corresponding 
to an integrated luminosity of 57.4 pb~^, taken in the years 1999-2000, at an ep centre of mass 
energy ^/s = 319 GeV, with a proton beam energy of 920 GeV. 

Only a short description of the HI detector is given here; a full description may be found 
in [25]. A right handed coordinate system is employed at HI that has its z-axis pointing in the 
proton beam, or forward, direction and x (y) pointing in the horizontal (vertical) direction. 
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Charged particles are measured in the central tracking detector (CTD). This device consists 
of two cylindrical drift chambers interspersed with ^-chambers to improve the ^-coordinate 
reconstruction and multi-wire proportional chambers mainly used for triggering. The CTD 
is situated in a uniform 1.15 T magnetic field, enabling momentum measurement of charged 
particles over the polar angular range' 20° < 9 < 160°. 

The CTD tracks are linked to hits in the vertex detector (central silicon tracker CST) [26] 
to provide precise spatial track reconstruction. The CST consists of two layers of double-sided 
silicon strip detectors surrounding the beam pipe, covering an angular range of 30° < 9 < 150° 
for tracks passing through both layers. The information on the ^-coordinate of the CST tracks 
is not used in the analysis presented in this paper. For CTD tracks with CST hits in both 
layers the transverse distance of closest approach (DCA) to the nominal vertex in x-y can be 
measured with a resolution of 33 fim © 90 yum/j»r[GeV], where the first term represents the 
intrinsic resolution (including alignment uncertainty) and the second term is the contribution 
from multiple scattering in the beam pipe and the CST; pt is the transverse momentum of the 
track. 

The track detectors are surrounded in the forward and central directions (4° < 6* < 155°) by 
a fine grained liquid argon calorimeter (LAr) and in the backward region (153° < 9 < 178°) by a 
lead-scintillating fibre calorimeter (SPACAL) [27] with electromagnetic and hadronic sections. 
These calorimeters provide energy and angular reconstruction for final state particles from the 
hadronic system. The SPACAL is used in this analysis to measure and identify the scattered 
positron. A planar drift chamber (BDC [28]), positioned in front of the SPACAL (151° < 9 < 
178°), measures the angle of the scattered positron and allows suppression of photoproduction 
background, where particles from the hadronic final state fake a positron signal. 

Electromagnetic calorimeters situated downstream in the positron beam direction allow de- 
tection of photons and electrons scattered at very low Q"^. The luminosity is measured from the 
rate of photons produced in the Bethe-Heitler process ep ep'j. 

4 Experimental Method 
4.1 Event and Track Selection 

The events are selected by requiring a compact electromagnetic cluster in the SPACAL associ- 
ated with a track segment in the BDC to define the scattered positron candidate. The z position 
of the interaction vertex, reconstructed by one or more charged tracks in the tracking detectors, 
must be within ±20 cm of the centre of the detector to match the acceptance of the CST. Pho- 
toproduction events are suppressed by requiring — j) > 35 GeV. Here, Ei and j 
denote the energy and longitudinal momentum components of a particle and the sum is over all 
final state particles including the scattered positron and the hadronic final state (HFS). The HPS 
particles are reconstructed using a combination of tracks and calorimeter deposits in an energy 
flow algorithm that avoids double counting. The event kinematics, and the inelasticity vari- 
able y, are reconstructed with the 'eS' method [29], which uses the scattered positron and the 

' The angular coverage of each detector component is given for the interaction vertex in its nominal position. 
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HFS. The Bjorken scaling variable x is obtained from x = j sy. In order to have good accep- 
tance in the SPACAL and to ensure that the HFS has a significant transverse momentum, events 
are selected in the range 6.3 < < 120 GeV^. The analysis is restricted to 0.07 < y < 0.7 
to ensure that the direction of the quark which is struck by the photon is mostly in the CST 
angular range. A further cut ofy< 0.63 is imposed for events with < 18 GeV^ to reduce 
photoproduction background. 

The triggers used in the analysis require a SPACAL energy deposit in association with a 
loose track requirement. Although these triggers are almost 100% efficient, not all events could 
be recorded, due to the large rate for low events. A fraction of events is rejected at the first 
trigger level (LI) and final trigger level (L4). The Monte Carlo events are assigned weights to 
account for the events rejected at LI while the data events are assigned weights to account for 
the events rejected at L4. The weights are largest for those events with an electron at low radius 
and low energy. The overall effect of the trigger is a reduction of the effective luminosity by a 
factor of about 10 for the lowest bin and 1.4 for the highest. After applying the event weights 
and the inclusive selection detailed above, the total number of events is about 1.5 million. The 
background from photoproduction events is estimated from the PHOJET Monte Carlo simula- 
tion. In most of the y range this background is negligible and does not exceed 9% in any x-Q'^ 
bin used in this analysis. 

The primary event vertex in r-0 is reconstructed from all tracks (with or without CST hits) 
and the position and spread of the beam interaction region [1]. The impact parameter of a track, 
which is the transverse distance of closest approach (DCA) of the track to the primary vertex 
point, is only determined for those tracks which are measured in the CTD and have at least two 
CST hits linked (referred to as CST tracks). Only CST tracks with a transverse momentum > 
0.5 GeV are included in the DCA and related distributions that are used to separate the different 
quark flavours. In the kinematic range of this measurement, the fraction of c (b) events that 
have at least one charged track within the angular range of the CST, with transverse momentum 
> 0.5 GeV and originating from the decay of a heavy flavoured hadron, is expected to be 82% 
(96%), as determined from the Monte Carlo simulation. The efficiency to obtain a CST track 
from a CTD track is 76%, within the angular range of the CST. 

In order to determine a signed impact parameter (5) for a track, the azimuthal angle of 
the struck quark 0quaik must be determined for each event. To do this, jets with a minimum 
Pt of 2.5 GeV, in the angular range 15° < 9 < 155°, are reconstructed using the invariant 

algorithm [30] in the laboratory frame using all reconstructed HFS particles. The angle 
0quark IS defined as the of the jet with the highest transverse momentum or, if there is no jet 
reconstructed in the event, as 180° — 0cicc, where 0cicc is the azimuthal angle of the electron in 
degrees. The direction defined in the transverse plane by 0quark and the primary vertex is called 
the quark axis. Approximately 81% (95%) of c (b) events have 0quark reconstructed from a jet, 
as determined from the Monte Carlo simulation. 

The difference between the reconstructed to the true (pqumk (defined as the azimuthal angle 
of the quark with highest transverse momentum) is estimated from the Monte Carlo simulation 
to have a resolution of about 5° for events with a reconstructed jet and 35° for the rest. The reso- 
lution of </)quark IS chcckcd with events containing a reconstructed D* meson. Figure[T]shows the 
difference between the reconstructed D* azimuthal angle and 0quark for events with and without 
a reconstructed jet. Both distributions are well described by the Monte Carlo simulation. 
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If the angle between the quark axis and the line joining the primary vertex to the point of 
DCA is less than 90°, 5 is defined as positive, and is defined as negative otherwise. Tracks with 
azimuthal angle outside ±90° of 0quark are rejected. The 6 distribution, shown in figure El is 
seen to be asymmetric with positive values in excess of negative values indicating the presence 
of long lived particles. It is found to be well described by the Monte Carlo simulation. Tracks 
with \5\ > 0.1 cm are rejected from the analysis to suppress light quark events containing long 
lived strange particles. 

4.2 Quark Flavour Separation 

The method used in [ 1] to distinguish between the c, b and light quark flavours has been modified 
in the present analysis because here the fraction of b quarks is smaller. The quantities Si, S2 
and 5*3 are defined as the significance {6 /a (6)) of the track with the highest, second highest and 
third highest absolute significance, respectively, where a (5) is the error on 5. Distributions of 
each of these quantities are made. The events contributing to the 5*2 distribution also contribute 
to the distribution. Similarly, those contributing to the S3 distribution also contribute to the 
S2 and distributions. Events in which Si and S2 have opposite signs are excluded from the 
5*2 distribution. Events in which Si, S2 and S3 do not all have the same sign are excluded from 
the 5*3 distribution. 

Figure |3] shows the three significance distributions. The simulation gives a reasonable de- 
scription of the data. In order to substantially reduce the uncertainty due to the resolution of 
5 and the light quark normalisation, the contents of the negative bins in the significance dis- 
tributions are subtracted from the contents of the corresponding positive bins. The subtracted 
distributions are shown in figure HI It can be seen that the resulting distributions are dominated 
by c quark events, with a b fraction increasing with significance. The light quarks contribute a 
small fraction for all values of significance. 

The fractions of c, b and light quarks of the data are extracted in each x-Q^ interval using a 
least squares simultaneous fit to the subtracted Si, S2 and S3 distributions (as in figure |4|) and 
the total number of inclusive events before any CST track selection. The c, b and uds Monte 
Carlo simulation samples are used as templates. The Monte Carlo c, b and uds contributions 
in each x-Q^ interval are scaled by factors Pc, Ph and Pi, respectively, to give the best fit to 
the observed subtracted 5*1, S2, S3 and total distributions. Only the statistical errors of the data 
and Monte Carlo simulation are considered in the fit. The fit to the subtracted significance 
distributions mainly constrains Pc and P^, whereas the overall normalisation constrains Pi. 

The results of the fit to the complete data sample are shown in figure IH The fit gives a 
good description of all the significance distributions, with a /n.d.f of 18.0/25. Values of 
Pc = 1.28±0.04, Pb = 1. 55 ± 0.16 and = 0.95 ±0.01 are obtained. The c and 6 scale factors 
are found to be anti-correlated with an overall correlation coefficient of -0.70. Acceptable 
values are also found for the fits to the samples in the separate x-Q"^ intervals. Since the same 
event may enter the 5*1, S2 and S3 distributions, it was checked using a high statistics Monte 
Carlo simulation that this has negligible effect on the results of the fits with the present data 
statistics. 
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The results of the fit in each x-Q'^ interval are converted to a measurement of the 'reduced 
c cross section' defined from the differential cross section as 



^'^^ dxdg2 27ra2(i + (1-^)2)' 

using: 

P N^'~^^^^ 

a'^'^fx, Q^) = Q^) rjr^ — ^77=; 07^ — ^bcc, (2) 

V !^ / V 7 / p ^MCgen _|_ p^jyMCgen _|_ p^jyMCgen ^ ^ 

where a{x, Q"^) is the measured inclusive reduced cross section from HI [28] and A*'^*-'^'^'^, 
are the number of c, h and light quark events generated from the Monte 
Carlo in each bin. A bin centre correction 5bcc is applied using a NLO QCD expectation for 
a'"^ to convert the bin averaged measurement into a measurement at a given x-Q"^ point. The 
NLO QCD expectation is calculated from the results of a fit similar to that performed in [3 1] but 
using the FFNS scheme to generate heavy flavours. A small correction (< 2.6%) for the beam 
energy difference is applied, using the NLO QCD expectation, to the measurement of cr(x, Q^) 
which was performed at a lower centre of mass energy of 301 GeV than the data presented here. 
The cross section is defined so as to include a correction for pure QED radiative effects. Events 
that contain c hadrons via the decay of h hadrons are not included in the definition of the c cross 
section. The differential h cross section is evaluated in the same manner. 



4.3 Systematic Errors 

The systematic uncertainties on the measured cross sections are estimated by applying the fol- 
lowing variations to the Monte Carlo simulation: 

• An uncertainty in the 5 resolution of the tracks is estimated by varying the resolution by 
an amount that encompasses the differences between the data and simulation (figures |2l 
Ol). This was achieved by applying an additional Gaussian smearing in the Monte Carlo 
of 200 /im to 5% of randomly selected tracks and 25 fim. to the rest. 

• A track efficiency uncertainty of 2% due to the CTD and of 1% due to the CST 

• The uncertainties on the various D and B meson lifetimes, decay branching fractions and 
mean charge multiplicities are estimated by varying the input values of the Monte Carlo 
simulation by the errors on the world average measurements. For the branching fractions 
of h quarks to hadrons and the lifetimes of the D and B mesons the central values and 
errors on the world averages are taken from [32]. For the branching fractions of c quarks 
to hadrons the values and uncertainties are taken from [33], which are consistent with 
measurements made in DIS at HERA [34]. For the mean charged track multiplicities the 
values and uncertainties for c and h quarks are taken from Marklll [35] and LEP/SLD [36] 
measurements, respectively. 

• An uncertainty on the fragmentation function of the heavy quarks is estimated using the 
Peterson fragmentation function [37] with parameters Cc = 0.058 and e?, = 0.0069, in- 
stead of the LUND fragmentation model. 
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• An uncertainty on the QCD model of heavy quark production is estimated by using the 
CASCADE Monte Carlo instead of the RAPGAP Monte Carlo. 

• The uncertainty on the asymmetry of the light quark 6 distribution is estimated by repeat- 
ing the fits with the subtracted light quark significance distributions (figurelH) changed by 
ib50%. The light quark asymmetry was checked to be within this uncertainty by compar- 
ing the asymmetry of Monte Carlo events to that of the data, in the region 0.1 < \6\ < 
0.5 cm, where the light quark asymmetry is enhanced. 

• An error on the quark axis is estimated by shifting the quark axis by 2°(5°) for events with 
(without) a reconstructed jet. These shifts were estimated by comparing the difference 
between 0quark and the track azimuthal angle in data and Monte Carlo. 

• A 4% uncertainty on the hadronic energy scale. 

• Uncertainties on the acceptance and bin centre correction due to the input structure func- 
tions used are estimated by reweighting the input a^'' distribution by and 1 ± 
0.21n[QV(10 GeV^)] and a^'^ by x^^-^ and 1 ± 0.41n[gV(10 GeV^)]. The range of 
variation of the input structure functions was estimated by comparing to the measured 
values obtained in this analysis. 

• An uncertainty on the photoproduction background is estimated by assigning ±100% of 
the expected number of events from the PHOJET simulation that enter the significance 
distributions. 

Other sources of systematic error pertaining to the NC selection were also considered [28]: a 
1.5% uncertainty on the luminosity measurement; an uncertainty on the scattered positron polar 
angle of 0.3 mrad and energy of 0.3-1.0% depending on the energy; a 0.5% uncertainty on the 
scattered positron identification efficiency; a 0.5-2% uncertainty on the positron track-cluster 
link efficiency; a < 1% uncertainty on the trigger efficiency and a 1% uncertainty on the cross 
section evaluation due to QED radiative corrections. 

A detailed list of the systematic effect on each cross section measurement is given in table[T] 
The systematic error is larger for the b measurement than it is for the c because the b fraction 
is much smaller than the c fraction. The errors which contribute most to the uncorrelated sys- 
tematic error in table [T] are, at low and high y, the uncertainty on the photoproduction 
background and, elsewhere, the uncertainty on the acceptance and bin centre correction due to 
the input structure function. 



5 Results 

The measurements of a'^'^ are listed in table [T] and shown in figure |5] as a function of x for fixed 
values of Q^. The HI data for a'^'^ are compared with the results extracted from D* meson 
measurements by HI [5] and ZEUS [6] obtained using a NLO program [38] based on DGLAP 
evolution to extrapolate the measurements outside the visible D* range. The measurements for 
a^'' from the present analysis and the D* extraction methods are in good agreement. 
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The a^'^ data are compared with two VFNS predictions from NLO QCD (see section im 
from MRST [2] and CTEQ [3], and with predictions based on CCFM [15] parton evolution. 
The predictions provide a reasonable description of the present data. 

The measurements of 5"^^ are also listed in table[l]and are shown in figure|^as a function of 
X for fixed values of Q^. This is the first measurement of ct^^ in this kinematic range. The a^'^ 
data are also compared with the two VFNS NLO QCD predictions and the CCFM prediction. 
The difference between the two VFNS NLO QCD calculations, which reaches a factor 2 at the 
lowest and x, arises from the different treatments of threshold effects by MRST and CTEQ. 
Within the current experimental errors these differences cannot be resolved. 

The structure function F2'^ is evaluated from the reduced cross section 



,.2 

~ cc 



a- = Fr--— (3) 



l + il-yf 

where the longitudinal structure function F^'^ is estimated from the same NLO QCD expectation 
as used for the bin centre correction. The structure function Fn'^ is evaluated in the same manner. 



The measurements F2^ and F2'' are shown as a function of in figure 17] and figure [HI The 
measurements of Fg^'^ and F2'' show positive scaling violations which increase with decreasing 
X. The data are compared with the VFNS QCD predictions from MRST and CTEQ at NLO 
and a recent calculation at NNLO [13]. The charm data are more precise than the spread in 
predictions of the QCD calculations. 

The measurements are also presented in table |2l and figure |9l in the form of the fractional 
contribution to the total ep cross section 

^2 cc j2 

r' = — - — /— — — . (4) 

The h fraction /^^ is defined in the same manner. In the present kinematic range the value of f^" 
is around 24% on average and increases slightly with increasing and decreasing x. The value 
of increases rapidly with from 0.4% at = 12 GeV^ to 1.5% at = 60 GeV^. The 
NLO QCD predictions of MRST shown in figure 13 are found to describe the data reasonably 
well. 



6 Conclusion 



The differential charm and beauty cross sections in Deep Inelastic Scattering are measured at 
low and Bjorken x using the impact parameters of tracks from decays of long lived c and 
h hadrons as reconstructed from the vertex detector. This is the first measurement of F^ in 
the low kinematic region. In this kinematic range the charm cross section contributes on 
average 24% of the inclusive ep cross section, and the beauty fraction increases from 0.4% at 
= 12 GeV^ to 1.5% at = 60 GeV^. The cross sections and derived structure functions 
F2^ and F2'' are found to be well described by predictions of perturbative QCD. 



11 



Acknowledgements 



We are grateful to the HERA machine group whose outstanding efforts have made this experi- 
ment possible. We thank the engineers and technicians for their work in constructing and main- 
taining the HI detector, our funding agencies for financial support, the DESY technical staff for 
continual assistance and the DESY directorate for support and for the hospitality which they ex- 
tend to the non-DESY members of the collaboration. We are grateful to S. Kretzer, R. S. Thorne 
and W. K. Tung for providing us with their calculations and for productive discussions. 



References 



[1] A. Aktas et al. [HI Collaboration], Eur. Phys. J. C 40 (2005) 349 | |hep-ex/0411046t 



[2] A. D. Martin, R. G. Roberts, W. J. Stirling and R. S. Thome, Eur. Phys. J. C 39 (2005) 155 
[ |hep-ph/0411040| . 

[3] S. Kretzer, H. L. Lai, E I. Olness and W. K. Tung, Phys. Rev. D 69 (2004) 114005 
[hep-ph/0307022). 

[4] C. MXoiietal. [HI Collaboration], Z. Phys. C 72 (1996) 593 | |hep-ex/9607(7n | ; 

J. Breitweg et al. [ZEUS Collaboration], Phys. Lett. B 407 (1997) 402 ||hep^ex/9706009J; 
C. AdloE etal. [HI Collaboration], Nucl. Phys. B 545 (1999) 21 [hep-ex/98 12023] ; 
J. Breitweg et al. [ZEUS Collaboration], Eur. Phys. J. C 12 (2000) 35 [ |hep-ex/99080T2| . 

[5] C. KdloE etal. [HI Collaboration], Phys. Lett. B 528 (2002) 199 | |hep-ex/0 1080391 . 



[6] S. Chekanov et al. [ZEUS Collaboration], Phys. Rev. D 69 (2004) 012004 
piep^ex/03080681. 



[7] S. Chekanov etal. [ZEUS Collaboration], Phys. Lett. B 599 (2004) 173 [ hep-ex/0405069 |. 

[8] A. Aktas et al. [HI Collaboration], Eur. Phys. J. C 41 (2005) 453 piep-ex/0502OTm . 

[9] E. Laenen, S. Riemersma, J. Smith and W. L. van Neerven, Nucl. Phys. B 392 (1993) 162; 
E. Laenen, S. Riemersma, J. Smith and W. L. van Neerven, Nucl. Phys. B 392 (1993) 229. 

[10] E I. Olness and W. K. Tung, Nucl. Phys. B 308 (1988) 813; 

M. A. G. Aivazis, E I. Olness and W. K. Tung, Phys. Rev. D 50 (1994) 3085 
[hep-ph/9312318|; 

M. A. G. Aivazis, J. C. Collins, E I. Olness and W. K. Tung, Phys. Rev. D 50 (1994) 3102 
[hep-ph/9312319|; 

M. Kramer, E I. Olness and D. E. Soper, Phys. Rev. D 62 (2000) 096007 
[he p-ph/0003035J . 



12 



[11] R. S. Thome and R. G. Roberts, Phys. Rev. D 57 (1998) 6871 | hep'-ph/9709442 1 ; 
R. S. Thorne and R. G. Roberts, Phys. Lett. B 421 (1998) 303 |hep-ph/971 1223|; 
R. S. Thorne and R. G. Roberts, Eur. Phys. J. C 19 (2001) 339 [hep-ph /0010344| ; 

[12] M. Buza, Y. Matiounine, J. Smith and W. L. van Neerven, Phys. Lett. B 411 (1997) 211 
[hep-ph/9707263|; 

A. Chuvakin, J. Smith and W. L. van Neerven, Phys. Rev. D 61 (2000) 096004 
[hep-ph/9910250|; 

A. Chuvakin, J. Smith and W. L. van Neerven, Phys. Rev. D 62 (2000) 036004 
[hep-ph/00020lT|. 

[13] R. Thome, "A Variable Flavour Number Scheme at NNLO," To appear in the proceed- 
ings of 13*'^ International Workshop on Deep Inelastic Scattering (DIS 2005), Madison, 
Wisconsin, USA, April 2005. [hep-ph/0506251 1. 

[14] M. Hansson and H. Jung, "Status of CCFM: Un-integrated gluon densities," Proceedings 
of ll*'^ Intemational Workshop on Deep Inelastic Scattering (DIS 2003), St. Petersburg, 
Russia, April 2003, p 488. Edited by V.T. Kim and L.N. Lipatov. |hep-ph/0309009|. 

[15] M. Ciafaloni, Nucl. Phys. B 296 (1988) 49; 

S. Catani, E Eiorani and G. Marchesini, Phys. Lett. B 234 (1990) 339; 
S. Catani, E Eiorani and G. Marchesini, Nucl. Phys. B 336 (1990) 18; 

G. Marchesini, Nucl. Phys. B 445 (1995) 49 | |hep-ph/9412327| . 

[16] H. Jung and G. P Salam, Eur. Phys. J. C 19 (2001) 351 |hep-ph/00l2l43|; 

H. Jung, Comput. Phys. Commun. 143 (2002) 100 [hep-ph/0109102|; 
(see also http://www.desy.de/~jung/cascade/(l. 

[17] H. Jung, Comput. Phys. Comm un. 86 (1995) 147; 
(see also |http ://www. desy.de/~ j ung/rapgap/ ) . 

[18] K. Charchula, G. A. Schuler and H. Spiesberger, Comput. Phys. Commun. 81 (1994) 381; 
(see also http://www.desy.de/~hspiesb/django6.html). 

[19] A. D. Martin, W. J. Stirling and R. G. Roberts, Phys. Rev. D 50 (1994) 6734 
[ hep-ph/9406315J . 

[20] M. Gliick, E. Reya and A. Vogt, Z. Phys. C 67 (1995) 433. 

[21] T. Sjostrand etal., Comput. Phys. Commun. 135 (2001) 238 [hep-ph/OOlOOT^. 

[22] A. Kwiatkowski, H. Spiesberger and H. J. Mohring, Comput. Phys. Commun. 69 (1992) 
155. 

[23] R. Engel and J. Ranft, Phys. Rev. D 54 (1996) 4244 | |hep-ph/9509373 |. 

[24] R. Brun, R. Hagelberg, M. Hansroul and J. C. Lassalle, CERN-DD-78-2-REV. 



13 



[25] 1. Abt et al. [HI Collaboration], Nucl. Instrum. Meth. A 386 (1997) 310 and 348. 

[26] D. Pitzl et al, Nucl. Instrum. Meth. A 454 (2000) 334 [ .hep-ex/0002044 |. 

[27] T. NichoUs et al. [HI SPACAL Group], Nucl. Instrum. Meth. A 374 (1996) 149. 

[28] C. Adloff a/. [HI Collaboration], Eur. Phys. J. C 21 (2001) 33 [hep-ex/00 12053}. 

[29] U. Bassler and G. Bemardi, Nucl. Instrum. Meth. A 361 (1995) 197 [hep-ex/94 12004} ; 
U. Bassler and G. Bemardi, Nucl. Instrum. Meth. A 426 (1999) 583 | |hep-ex/98010T7| . 

[30] S. D. Ellis and D. E. Soper, Phys. Rev. D 48 (1993) 3160 |hep-ph/9305266}; 

S. Catani, Y. L. Dokshitzer, M. H. Seymour and B. R. Webber, Nucl. Phys. B 406 (1993) 
187. 

[31] C. Kdloiietal. [HI Collaboration], Eur. Phys. J. C 13 (2000) 609 [hep-ex/9 908059 |. 

[32] K. Hagiwara et al. [Particle Data Group], Phys. Rev. D 66 (2002) 010001. 

[33] L. Gladilin, "Charm hadron production fractions", piep^x/9912064|. 

[34] A. Aktas et al. [HI Collaboration], Eur. Phys. J. C 38 (2005) 447 [hep-ex/0 408149| . 

[35] D. Coffman et al. [MARK-III Collaboration], Phys. Lett. B 263 (1991) 135. 

[36] D. Abbaneo et al. [LEP/SLD Heavy Flavour Working Group], LEPHF 2001-01; 
(available from http://lepewwg.web.cern.ch/LEPEWWG/heavy/). 

[37] C. Peterson, D. Schlatter, I. Schmitt and PM. Zerwas, Phys. Rev. D 27 (1983) 105. 

[38] B. W. Harris and J. Smith, Phys. Rev. D 57 (1998) 2806 | |hep-ph/9706334| |. 



14 







X 


y 




Cch 


(^stat 


^sys 


5tot 








<^Dmul 




'^frag 


'^modcl 


^uds 




Ff 




(GeV2) 


(•10-3) 








(%) 


(%) 


(%) 


(%) 


(%) 


(%) 


(%) 


(%) 


(%) 


(%) 


(%) 


(%) 




c 


12 


0.197 


0.600 


0.412 


-0.62 


12 


13 


18 


11 


+3.2 


-1.4 


-3.1 


-0.3 


-0.7 


-1.9 


-5.0 


+2.0 


0.435 


c 


12 


0.800 


0.148 


0.185 


-0.68 


8.8 


9.4 


13 


5.6 


+2.5 


-1.7 


-3.2 


-0.2 


-0.4 


-2.2 


-5.2 


+2.0 


0.186 


c 


25 


0.500 


0.492 


0.318 


-0.66 


8.7 


10 


13 


6.8 


+3.1 


-1.4 


-3.1 


-0.3 


-0.7 


-1.9 


-5.0 


+2.0 


0.331 


c 


25 


2.000 


0.123 


0.212 


-0.72 


5.2 


8.6 


10 


4.1 


+2.6 


-1.6 


-3.1 


-0.2 


-0.5 


-2.1 


-5.2 


+2.0 


0.212 


c 


60 


2.000 


0.295 


0.364 


-0.74 


6.2 


8.3 


10 


3.5 


+3.2 


-1.4 


-3.1 


-0.3 


-0.7 


-1.9 


-5.0 


+2.0 


0.369 


c 


60 


5.000 


0.118 


0.200 


-0.76 


7.8 


8.5 


12 


3.8 


+2.7 


-1.6 


-3.1 


-0.2 


-0.5 


-2.1 


-5.1 


+2.0 


0.201 


& 


12 


0.197 


0.600 


0.0045 


-0.62 


55 


22 


60 


12 


-13 


-7.5 


-2.9 


+3.0 


+4.6 


+8.9 


-4.8 


+1.3 


0.0045 




12 


0.800 


0.148 


0.0048 


-0.68 


30 


33 


45 


13 


-21 


-10 


-5.4 


+3.1 


+6.9 


+15 


-7.7 


+1.7 


0.0048 




25 


0.500 


0.492 


0.0122 


-0.66 


22 


21 


31 


9.1 


-13 


-7.6 


-3.0 


+3.0 


+4.7 


+9.1 


-4.8 


+1.3 


0.0123 


& 


25 


2.000 


0.123 


0.0061 


-0.72 


26 


28 


39 


9.8 


-18 


-9.4 


-4.7 


+3.1 


+6.3 


+13 


-6.8 


+1.6 


0.0061 


b 


60 


2.000 


0.295 


0.0189 


-0.74 


21 


20 


29 


6.2 


-13 


-7.5 


-2.9 


+3.0 


+4.6 


+8.8 


-4.7 


+1.3 


0.0190 


b 


60 


5.000 


0.118 


0.0130 


-0.76 


26 


25 


36 


7.4 


-16 


-8.8 


-4.1 


+3.0 


+5.8 


+12 


-6.1 


+1.5 


0.0130 



Table 1: The measured reduced NC cross section (a'^^) for charm (c) and beauty (b) quarks, shown with the correlation coefficients 
(Ccb), the statistical error i5sta.t), the systematic error {5sys), the total error ((5tot) and the uncorrelated systematic error (5unc)^ The next 8 
columns represent a +lcr shift for the correlated systematic error contributions from: track impact parameter resolution, track efficiency, 
D multiplicity, B multiplicity, fragmentation, QCD model, light quark contribution and quark axis 0quark^ The —la errors are taken as 
the negative of the upward errors. The errors are correlated between charm and beauty but uncorrelated to inclusive data, apart from a 
normalisation uncertainty of 1.5% which is 100% correlated. The table also shows the values for Fg'^ and F2'' obtained from the measured 
cross sections using the NLO QCD fit to correct for the contributions from F^^ and The quoted relative errors apply also to Fg'^ and 
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Table 2: The measured charm (f'^'^) and beauty (/ ) fractional contributions to the total ep cross 
section, shown with statistical (^staf ^stat)' systematic (6^ys, S^L) and total (6^^^^, 6^^^) errors. 




Figure 1: The azimuthal difference between the D* and the quark axis for those events where 
the quark axis is defined (a) by a jet and (b) by 180° — 0cicc- Included in the figure is the 
expectation from the Monte Carlo simulation normalized to the number of data events. 
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Figure 2: The distribution of the signed impact parameter 5 of a track to the primary vertex in the 
x-y plane. Included in the figure is the expectation from the DJANGO Monte Carlo simulation 
for light quarks and that from the RAPGAP Monte Carlo simulation for c and h quarks. The 
contributions from the various quark flavours are shown after applying the scale factors obtained 
from the fit to the subtracted significance distributions of the data (see section lOl) . 
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Figure 3: The significance 5/cr((5) distribution (a) of the highest absolute significance track 
(^i), (b) of the track with the second highest absolute significance (52) and (c) of the track with 
the third highest absolute significance {S3). Included in the figure is the expectation from the 
DJANGO Monte Carlo simulation for light quarks and that from the RAPGAP Monte Carlo 
simulation for c and b quarks. The contributions from the various quark flavours are shown 
after applying the scale factors obtained from the fit to the subtracted significance distributions 
of the data. 
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Figure 4: The subtracted significance distributions of (a) 5*1, (b) 5*2 (c) S3. Included in the 
figure is the result from the fit to the data of the Monte Carlo distributions of the various quark 
flavours. 
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Figure 5: The measured reduced cross section a^^ shown as a function of x for 5 different 
values. The inner error bars show the statistical error, the outer error bars represent the statistical 
and systematic errors added in quadrature. The measurements of a^'^ from HI at high values of 
[1], the measurements obtained from D* mesons from HI and ZEUS [5, 6] and predictions 
of QCD are also shown. 20 




Figure 6: The measured reduced cross section a^^ shown as a function of x for 5 different 
values. The inner error bars show the statistical error, the outer error bars represent the statistical 
and systematic errors added in quadrature. The measurements of a^^ from HI at high values of 
[1] and predictions of QCD are also shown. 
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Figure 7: The measured shown as a function of Q'^ for various x values. The inner error 
bars show the statistical error, the outer error bars represent the statistical and systematic errors 
added in quadrature. The F2'^ from HI at high values of Q"^ [1], the measurements obtained 
from D* mesons from HI and ZEUS [5, 6] and predictions of QCD are also shown. 
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Figure 8: The measured F!^'^ shown as a function of for various x values. The inner error 
bars show the statistical error, the outer error bars represent the statistical and systematic errors 
added in quadrature. The from HI at high values of [1] and predictions of QCD are also 
shown. 
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Figure 9: The contributions to the total cross section f^'^ and /^^ shown as a function of for 
6 different x values. The inner error bars show the statistical error, the outer error bars represent 
the statistical and systematic errors added in quadrature. A prediction of NLO QCD is also 
shown. 
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